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■ Non-perturbative QCD vacuum effects at finite density are parametrized by means of a bag con- 

^Nj ' stant B. It is extracted fi'om a Nambu- Jona Lasinio model with two or three flavors. The parameter 

B is used in an effective bag-hke model of baryons to study the nuclear phenomenology. We examine 
(N . the nucleon structure and the thermodynamical properties of symmetric nuclear matter, particular 

I attention is paid to the symmetry energy and to the eventual phase transition to deconfined quark 

matter. An alternative sketch of the binding mechanism of symmetric nuclear matter emerges within 
this approach. It is also found that the inclusion of strangeness content in B is crucial for an appro- 
priate description. 
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Abstract 



PACS:12.39.Ba,12.39.Fe,21.65.+f, 12.40.Yx, 21.30.Fe 



The low energy regime of QCD is very difficult to manage due to its non-perturbative nature. This 
\ feature gives rise to different methods, like sum rules and effective models, to avoid dealing with the 

full theory in the description of the hadronic phenomenology. These approaches intend to include the 
symmetries of QCD implicitly into the final results. 

A simple shortcut procedure to take care of non-perturbative vacuum contributions is the introduction 
of the so-called bag constant B, in both confined and deconfined quark phases. The former bag models 
used it to account for the energy required to create a unit volume of the Weyl-Wigner vacuum into the 
complex QCD ground state. In ref.^ two different scales were associated to B, the higher one is related 
to the condensation of gluons and consequently with the recovery of the scale symmetry of QCD. The 
lower one is indicative of the restoration of chiral symmetry. Approximate calculations for these scales 
there given are S^/"* = 24:5MeV and B^^^ — liOMeV respectively. 

Another significant aspect in the definition of the bag parameter is its hypothetical variation with 
the medium properties of the hadronic matter. Speculation about the medium dependence of B has 
been motivated from different areas and ad-hoc phenomenological parametrizations have been proposed 
[3 El EI El El • It is expected that the gradual recovery of fundamental symmetries with increasing density 
and/or temperature modifies the vacuum structure and the value of B accordingly. 

In this work we attempt to extract a bag parameter that varies with the baryonic density but carrying 
the signature of the chiral symmetry and, ultimately to examine its infiuence over the nuclear dynamics 
at finite density. For this purpose we use two rather different theoretical descriptions of the high density 
regime of the strong interacting matter and we combine them to describe in a simplified manner the 
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influence of the QCD vacuum changes over the gross properties of the nuclear matter. We have selected 
the Nambu-Jona Lasinio (NJL) model [7] to study the modification of the QCD vacuum with increasing 
baryonic density and to extract a medium-dependent bag constant. This parameter is then inserted 
into the Quark Meson Coupling (QMC) model [H |^ to describe the nucleon substructure. It must be 
mentioned that each one of these issues have been independently investigated in the past, but we intend 
to put them together in a novel combination. For instance, a density dependent bag constant has been 
proposed within the QMC 2 , but it was initially regarded as a way to restore certain aspects of the 
nuclear phenomenology instead of a manifestation of the underlying QCD vacuum. This fact explains 
why the parametrization of B has been made in terms of either the scalar a meson mean field or the 
effective nucleon mass. 

Adopting the viewpoint of the early bag models, we expect that B would carry significative information 
about the fundamental symmetries and would transfer it to the hadronic sector. Thus, as previously 
mentioned, the recovery of chiral and scale symmetries with increasing density and/or temperature would 
have remarkable effects over B. 

The formulation of chiral invariant theories was certainly one of the major successes in the development 
of the bag models [23 • However QMC abandoned this requirement since it was thought as a way to match 
the bag picture of the hadron structure with the Quantum Hadrodynamics treatment of the nuclear 
many-body problem. The coupling to a scalar a meson plays a key role in this description, therefore the 
explicit breaking of chiral symmetry seems to be an inherent feature of the QMC formulation. In the 
following we shall see how a bag constant respecting the chiral invariance in its definition, could modify 
this situation. 

Our starting point is the original meaning of B, i.e. the energy difference between the fundamental 
state with explicit realization of the symmetries and that with symmetries spontaneously broken ^21 ■ For 
this purpose we consider the NJL for both SU(2) and SU(3) versions of the chiral symmetries |13[I14[[T5] . 
In the last case the fiavor degeneracy is removed introducing a heavy strange quark. 

The lagrangian density of the model is JS] 

C^q{ip-m)q + Ci + Ce (1) 

with g* = {u,d,s) and m =diag(mi,, rritj, ms) is the current mass matrix that explicitly breaks chiral 
invariance. The four quarks interaction term £4 is invariant under Ul{3) x Ur{3) transformations 

£4 - G [{qXkqf + iqil5Xkqf] , (2) 

where Xk, k — 0, 8 are the Gell-Mann matrices with Aq — ■\/2/3 diag(l, 1, 1). 
The six quarks term removes the redundant C/a(1) symmetry 

£e = -i^{det[g,(l-H75)gj] +det[g,(l-75)qj]} (3) 

Due to the development of the quark-antiquark condensates the current quark masses acquires its 
large constituent values Mk, given by the self-consistent equations jl5| 

Mk ^mk-4:G < qtqk > +2K < q,q, >< q^q^ >, (4) 

with i ^ k =/= j , and the condensates are given by 

<qkqk> = Ck^-i lim f ^ e'^P-^^' TrGkip) (5) 



where we have used the quark propagator 
Gkip) = ^+Mk) 



AEa - PkEp, - Ml ln( A±|^) 

Pk + tL„ 



p2 - 



^ - + 2m5{p'-Ml)e{pk~p) 



(6) 



(7) 



and Ex = yjx^ + M'^. We have used a non-covariant 3-momentum cutoff A, and explicitly included finite 
density effects using the Fermi momentum pk- Therefore the ground state is obtained filling the quark 
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Fermi shell plus the Dirac sea divergent contribution. The contribution of fe-flavor to the total baryon 
density is related to its Fermi momentum by 



n. = ^||. (8) 

Once the eqs. Q and © have been solved the energy density for the fundamental state is evaluated 
by taking the statistical average of the hamiltonian density 

^Q^' = ^ Y.^F{Mk,Pk) - F{Mk, A)) + 2G{Cl + Cl + Cl) ~ AKCuCaCs, (9) 

k 

with F{M,x) = xE^{El + - M^hi{{x + E^)/M). 

A definite constant Eq must be added to Eq. Q in order to get zero energy when no valence quarks 
are present. The pressure is evaluated using the thermodynamical relation 

Pqm = ^ ^^knk - eqm, (10) 
fc 

where /Ltfc = Osqm /duk is the chemical potential for quarks of flavor k. 

As above mentioned we take B as the difference between Eq. and the energy of the state with 
zero quark condensates, but at the same baryonic density. This definition is in agreement with that 
used in ^H]- In this way a density dependent bag parameter has been obtained, able to be used in both 
hadron and quark matter since the parametrization is made in terms of only a common conserved charge. 
It is expected that the flavor composition of matter affects the density dependence of B, therefore we 
examine three different situations (a) symmetric u — d quark matter in the SU(2) version of the model, (b) 
symmetric u — d quark matter in the SU(3) version, that is the Dirac sea contribution of strange quarks is 
included, and (c) electrically neutral quark matter coexisting with leptons in equilibrium by electroweak 
decay, a situation of interest in astrophysical applications. Thus we have the constraints for the quark 
Fermi momenta pu = Pd = P,Ps = 0, n = 2NcP^ / (97r^) for the items (a)-(b) and Hs = l^d, fJ'u + IJ'e = fJ-d, 
= Nc{2pl^ ~ Pd" pI)I^~ pI^ ^c{pt+Pd+ PD/i^'^'^) for the case (c). Here n stands for the baryon 
number density and He for the electron chemical potential. 

For practical calculations we have used = rud = 5.5MeV, — lA0.7AIeV and the model param- 
eters A = 602.3 MeV, GA^ = 1.835, KA^ = UMMeV which were obtained by reproducing the tt-, K- 
and Ty'-meson masses together with the constant ^31) for cases (b) and (c). For the item (a) instead, 
we have used the parameters given in |15) for the Hartree results in the 3-momentum cutoff approach, 
i.e. A = 653MeV, GA^ = 2.14. 

The density dependence obtained in our calculations is shown in Fig.l. The SU(2) and SU(3) schemes 
differ essentially in all the range, while within the three flavor treatment the cases (b) and (c) are similar 
for low and medium densities. The items (a) and (c) exhibit a steep descent at sufficiently high densities, 
in contrast with case (b) that keeps a large constant value. The cause of this behavior is that the strange 
condensate does not melt in (b) as the non-strange ones do. In turn this is due to the lack of the valence 
strange quark contribution, as confirmed by the sudden splitting of the curves (b) and (c) at n ~ 4.5 no 
shown in Fig. 1. The fast decrease of the graph corresponding to (c) coincides with the turning up of the 
strange quark into the Fermi sea. Therefore the large asymptotic value shown by the SU(3) approach for 
the symmetric u — d matter can be assigned to the constraint of null strange Fermi momentum. 

The resulting B parameter can be used to evaluate nuclear matter properties. For this purpose we have 
selected a bag-like model that provides a satisfactory description of a wide range of nuclear phenomena, 
the QMC model 1^1 ■ Within the QMC baryons are regarded as non-overlapping spherical bags, where 
three valence quarks are confined. These quarks interact with its surrounding media by the interchange 
of light (T and w-mesons, as suggested by the successful Quantum Hadrodynamics models |11| . Inside 
the bag the quark fields obey the mean-field equation 

ii^d^ ~ gtl'coo - mk*)qk{x) = 0, (11) 

where = — g^a, since within the bag there is only a small breakdown of chiral symmetry we use 
the current value of the quark masses . In Eq. pi|l a and loq stand for the mean field values of the 
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meson fields evaluated in the dense hadronic medium and 5, 
The solutions of this equation is |H] 



fc 



are the quark-meson coupling constants. 



An 



jo{xk r/R) 
i(3kS-fji{xk r/R) 



(12) 



where x is the quark spinor and 

ek 



R 



k , k 



N = R^2n^{nk-l) + Rmk*] 



Joixk) 



f3k 



flh - Rml 



1/2 



Qk + Rml^ 

with Qk — [x^. + (i?m^)^]^/^. The eigenvalue Xk is solution of the equation 

joixk) = Pk ji{xk) 

which arises from the condition of zero quark current through the bag surface. 
In this model the ground state bag energy is identified with the baryon mass M* 



M* 



R 



4 , 
-.BR^ 



(13) 



(14) 



(15) 



where Nk is the number of quarks of flavor k inside the bag. As usual in the bag model, two phe- 
nomenological corrective terms are added to the quark energy: the one containing zq amends the single 
particle spectrum and the other, proportional to the bag volume, represents the energy required to create 
a Weyl-Wigner bubble in the complex QCD vacuum. The last one has acquired a density dependence in 
our approach. 

The bag radius adjusts to maintain the confinement volume in equilibrium with it surrounding media. 
Thus once B is given and zq fixed to reproduce the baryon mass at zero baryon density, the bag radius 
R is determined by the equilibrium condition dM* /dR = 

Density dependent bag parameters have been proposed previously within the QMC description |21E1^- 
Formerly it was regarded as a way to restore the phenomenology oi a — u) fields at the saturation nuclear 
density uq = 0.15/m~^ However the density dependence of B was guessed by means of the effective 
nucleon mass or by the cr-mean field, thus its relation with the underlying theory was not clear. Instead 
a phenomenological relation between B and the equation of state of nuclear matter was the proposal in 
[H], and restriction to experimental electron scattering analysis was imposed in The decrease of B at 
0.7no is of 55% and 22% for curves (a) and (b)-(c) respectively. In the last case it is only a few percent 
above the upper bound of 10 — 17% established in 0]. Based on different grounds, the necessity of a 
variable bag parameter was claimed in [S| in order to describe appropriately the structure of neutron 
stars. 

Next we shall examine the behavior of symmetric nuclear matter at medium and high densities under 
the varying B-parameter that includes the signatures of the chiral symmetry of the strong interaction, 
considering two or three fiavors contributions. 

The thermodynamical properties in the QMC model are evaluated regarding nucleons as quasi- 
particles with single particle spectrum = ^/p^ + M*^ + gujUJo, where the effective nucleon mass must 
be evaluated from Eq. ()15|l . The energy density and the pressure for symmetric nuclear matter are, 
therefore [H] 



Pnm = n€pp — Enm- 



(16) 
(17) 



In Eq. (UHl the energy of mesons has been included, p_F is the nucleon Fermi momentum related with the 
baryon density n = 2p'p/(37r^), and conditions of homogeneity and static equilibrium has been assumed. 
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The symmetry energy of nuclear matter is a very significative quantity to understand a variety of nuclear 
processes at medium and high densities 17 ^. For example, it has a crucial role in the dynamics and 
structure of neutron stars, as well as in the composition of matter in heavy-ion collisions. However its 
behavior far from the normal nuclear density is poorly known experimentally, and theoretical predictions 
range over a wide and sometimes contradictory qualitative results. For a linear nuclcon-p meson coupling 
the energy symmetry for symmetric nuclear matter is 



The mesons mean field values have not been determined yet, they can be obtained minimizing the 
energy density with respect to a and uj. This leads to the equations |5] 



The nuclcon- meson couplings are expressed in terms of the quark- meson couplings as = 3gJ, g^j — Sg^. 
Numerical values can be assigned reproducing the conditions of nuclear matter at the saturation density 
no- Meaningfully it is found in the SU(3) treatment that, for a certain value of the binding energy, one 
can do calculations without the ct- meson coupling, i.e. = 0. Thus we choose the values PNM{no) = 0, 
£NM{nQ) — Mq = — 16.807MeT^ for the pressure and the binding energy at the normal density. If the 
vector-isovector meson p is also considered, then the corresponding coupling constant gp = g^ can be 
deduced equating the symmetry energy esym — (deN/dn) at ng to its phenomenological value esym = 
34:MeV. The coupling constants obtained under these assumptions are shown in Table I. 

Numerical calculations have been developed taking — 2.787/m~^, to^^ — and irip = 

3.902 fm^^ for the meson masses, Alp = 4.755/r7i^^, A/„ = 4.761 fm^^ for the nucleon masses at zero 
baryon density. For the sake of comparison we have also carried out calculations with the standard fixed 
B treatment (in the following labelled case (d)). 

Some results for the in-medium nucleon properties are displayed in Fig. 2. The neutron mass is 
uniformly decreasing, the falling at no is approximately 10% in the SU(3) treatment, 20 — 25% in the 
constant B option and of an excessively high 40% in the SU(2) case. The behavior of the bag radius in 
terms of the baryon density is shown in Fig. 2, the cases (b) and (d) provide a radius with a bounded 
variation. In particular for the approach (b) an increasing radius is obtained whose growth saturates 15% 
over its zero density value. Instead for case (a) a monotonously increasing radius is obtained, exceeding 
the limit of validity of the model. The large increment of the bag radius yields a violation of the non- 
overlapping bag hypothesis at n ~ no, thus we may conclude that the two fiavor description of the bag 
parameter is not able of true physical consideration. In the three fiavor case this assumption is also 
broken down but at the higher density n ~ 6 no . 

The results for the pressure and the symmetry energy are shown in Fig. 3, it can be seen that the 
softest equation of state at high densities corresponds to (d) . The nuclear compressibility k = 9 BPnm / dn 
takes on the values 310 and 335Mey//m'^ at n — uq for the (d) and (b) approaches, respectively. For 
the symmetry energy all the results are very similar for low densities and differences becomes sensible 
only beyond n/no ~ 2.5. 

Thermodynamical quantities are similar in the approaches (b) and (d) , corresponding greater pressures 
and lower symmetry energy to the instance (b). Instead nucleon mass and radius have distinguishing 
behaviors in all the range of densities. 

The medium dependent B description has an upper bound of applicability given by the point where 
bags start to overlap, unless new physical aspects arise before this situation is reached. This could be 
the case if a deconfinement phase transition take place. To examine this possibility we have compared 
the equation of state obtained previously with that corresponding to quark matter. The Gibbs criterium 
for phase transitions Pnm = Pqm a-nd /i„ = 2fij^ + fi^ is adopted, and variable or fixed B is used. Thus 
a first order phase transition is found only for the case (a), but at so high density that QMC model 
assumptions have been violated. Not transition at all is found for the cases (b). 




(18) 




(19) 



(20) 
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In this paper we have studied the effects of including the chiral symmetry of QCD and the strange 
degree of freedom in the definition of the bag parameter B. This has been obtained from a Nambu- Jona 
Lasinio model with two or three flavors, at finite baryon density. Since baryon number is a conserved 
charge of both quark and hadron matter, the quantity B obtained can be regarded as medium-dependent 
and applied to examine nucleon properties and the nuclear matter equation of state. Thus our results for 
B are based on a dynamical effective model of QCD, in contrast with previous calculations that used ad 
hoc parametrizations in terms of the cr-meson mean field. The description in terms of only non-strange 
quarks must be discarded, i.e. the bag parameter must include the strange vacuum contribution even for 
the description of ordinary nuclear matter. Within the SU(3) scheme we have found that the behavior 
of B strongly depends on the composition of matter at medium and high densities. 

For the hadronic phase we used the QMC model of confined quarks coupled to mesons, that allows 
us to extract the confinement volume medium-dependence. From the comparison of the items (b) and 
(d) we conclude that the case (b) does not need of the a meson to appropriately describe the properties 
of nuclear matter. Therefore the saturation mechanism differs from those of Quantum Hadrodynamics, 
the attractive contribution is provided by the interplay between many-body effects and the modification 
of the hadronic phase vacuum. Avoiding a direct coupling between quarks and the cr-meson opens new 
possibilities such as a chiral extension of the QMC, or could eliminate the unpleasant features reported 
in [TBI. 

We have not found any first order phase transition from symmetric nuclear matter to quark matter, 
but this could be a failure of the approach that could be mended including quark correlations between 
neighboring bags in the high density regime of the hadronic phase |19| . 
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Case 


9a 


9lo 


9p 


Bl'^[MeV] 


i?o[/m] 


(a) 


11.97 


16.46 


8.43 


165.06 


0.83 


(b) 


0.00 


6.00 


9.53 


217.59 


0.57 


(dl) 


17.39 


8.26 


9.17 


165.06 


0.83 


(d2) 


18.30 


9.23 


9.04 


217.59 


0.57 



Tabic 1: Coupling constants for the QMC model for the cases considered: (a) and (b) correspond to 
medium dependent B obtained from a Nambu-Jona Lasinio model with SU(2) and SU(3) chiral symme- 
tries respectively, (dl) and (d2) correspond to the constant B fixed at the zero density value of (a) and 
(b) respectively. In the last two columns the zero density value of the bag parameter and the bag radius 
are shown. 
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Figure 1: The bag parameter in terms of the normahzed baryon number density. 
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Figure 2; In- medium nucleon properties as functions of the normalized baryon number density. In the 
upper figure the effective neutron mass relative to its value at zero baryon number density, in the lower 
figure the equilibrium bag radius. The line convention is explained in each figure, the labels a,b, and d 
have the same meaning as in the text. In the last case we have taken two different constant values for B 
obtained at zero baryon density in the treatments (a) and (b). 
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Figure 3: Thermodynamical properties of symmetric nuclear matter as functions of the normalized baryon 
number density. In the upper figure the symmetry energy and the pressure in the lower panel. For the 
line convention see the explanation in Fig. 2. 
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